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We investigate the scaling behavior of the B meson decay constant fs and fs s at (5=5.7, 5.9, 6.1, employing 
the NRQCD heavy quark action and the clover light quark action. Mixing effect from dimension-4 operator in the 
heavy-light axial-vector current is studied, and we find that the a dependence of fs is significantly reduced. Our 
preliminary result for the decay constants in the quenched approximation is /b=162(1 18 ) MeV, /b, = 190(1*9) 
MeV, and /B.//fl=1.18(lg). 



1. Introduction 

A recent development in the NRQCD study of 
heavy quarks on the lattice is the realization that 
the mixing of a dimension-4 operator with the 
axial- vector current, while nominally 0(a s a), has 
a significant effect in the value of the heavy-light 
decay constant p],0]. An investigation of how this 
mixing effect affects the scaling behavior of the 
decay constant is an important issue. 

In this work we study this problem, through 
simulations, with and without the operator mix- 
ing taken into account, at three values of (3. A 
comparison is also made of the present NRQCD 
results with our previous calculation with the rel- 
ativistic heavy quark action^. 

2. Method 

We describe the light quark by the O(a)- 
improved SW clover action with one-loop cor- 
rected c BW as in Ref. ||. For heavy quark, we 
employ two types of the NRQCD action and op- 



erator, one including all terms up to 0(1/M) and 
the other up to 0(l/M 2 ). 

The 0(1/M) NRQCD action we use is 



S = ^Q(t,x) 
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where Q is a two-component heavy quark field, 
H =-A (2) /[2M } and 6H=-gcr-B/[2M ]. To 
the same order in 1/M, the four-component Dirac 
field iph is related to Q via FWT transformation, 

(?<*>)■ (2) 

The mixing relation between the continuum 
axial-vector current and lattice counterparts, con- 
sistently expanded to 0(a s a) and 0(a s /M), is 
given by 



Aa 



l+aJ^}j^+a sP ^+a sP ^ 2 \ (3) 



* Presented by K-I. Ishikawa. 



where J^—'tpi^tph with tpi the light quark field 
and r=7574, J {1) =-^iY~faA. {±) ij; h and J (2) =^; 
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Table 1 

Parameters of simulation. 



/3 


6.1 


5.9 


5.7 


Vol. 


24 a x 64 


16 3 x 48 


12 3 x 32 


# of conf. 


120 


300 


300 


a" 1 [GeV] 


2.29 


1.60 


1.04 



7'aAf ± 'r^/ l . An important point observed in 
the first calculation of the one-loop coefficients 
/ 9 a' 1 ' 2 ^0 i s that the coefficient is not sup- 
pressed by 1/aM and remains as 0(1) for heavy 
quark, so that the mixing with the operator 
yields a large 0(a s aA QCD ) contribution. We have 
calculated the mixing coefficients for our 0(1/ M) 
NRQCD action which is slightly different from 
that of Ref. §. 

3. Results on mixing effects 

We carry out simulation at three values of /3 
employing lattices and statistics as listed in Table 
|l]. To set the lattice scale, we interpolate string 
tension data collected in Ref. [Q and set v / cx=427 
MeV. 

Figure |l| shows our results for the quantity <!> p 
= (a s (M P ) / a s (M B )) 2/ Po fp^/Mp~ at /3=5.9. We 
observe that the contribution of the mixing op- 
erators (0(a s a)), which is the difference between 
(O's) and (#'s) in the figure, is as large as that of 
the multiplicative renormalization of the leading 
operator (0(a s )), which is the difference between 
(O's) and (O's). This effect becomes more signif- 
icant towards heavier quark mass due to a large 

(2) 

value of p A and that of the matrix element of 
j( 2 \ so that the slope of <i>p becomes reduced 
with the inclusion of the mixing, as observed in 
Ref. |2|. We find this behavior to be more pro- 
nounced at /3 = 5.7. 

In Fig. H we compare results for 0(1/M) 
NRQCD action with previous JLQCD results 
obtained with the SW clover action for heavy 
quark, interpreted as a non-relativistic action 
within the Fermilab formalism ||, at /3=6.1. Since 
the latter calculation does not include the effect 
of 0(a s a) mixing, we plot NRQCD results for 
the one-loop corrected leading operator. A good 
agreement of results for the two actions provides 
a check of viability of both the 1/M expansion 
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Figure 1. 1/Mp dependence of $p with the 
0(1/M) NRQCD action at j3 = 5.9. 
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Figure 2. <J> P with the 0(1/ M) NRQCD action 
(this work) and with the SW clover action for 
heavy quark || at (3 — 6.1. Mixing effects are not 
included for both results. 

approach of NRQCD and the Fermilab interpre- 
tation of the clover action for heavy quark. 

Figure || presents the scaling behavior of Jb 
without (open symbols) and with (filled sym- 
boles) operator mixing, and for two choices of the 
momentum scale q*=7r/a and 1/a for the cou- 
pling constant. A large scatter of the values at 
a _1 wl GeV -1 (f3— 5.7) shows that one-loop es- 
timates of renormalization factors are not reli- 
able at such a large lattice spacing. This prob- 
lem is substantially alleviated at /3=5.9 and 6.1 
(0-4<a<0.6 GeV- 1 ). In this region, the NRQCD 
result without including the operator mixing con- 
tribution has a large a dependence, which is siz- 
ably reduced with full inclusion of the mixing. 

It is gratifying that the value of Jb in this range 
of a are reasonably consistent with the results 
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Figure 3. Scaling behavior of fs with the 0(1/M) 
NRQCD action. Results with the SW clover ac- 
tion for heavy quark || are also plotted. 

from the clover quark action (A's)|| over the 
same range. Strictly speaking, such a compari- 
son is to be made with the continuum extrapo- 
lated value of the latter. A mild scaling violation 
exhibited by the clover result suggests that the 
agreement would not be severely violated in such 
an extrapolation. Two points, however, have to 
be checked to consolidate the agreement: (i) the 
NRQCD values suffer from 0(a s /(aM) 2 ) errors 
toward smaller lattice spacing, whose magnitude 
in our range of a needs to be examined, (ii) The 
clover result does not incorporate effects of the 
0(a s a) mixing, whose magnitude is yet unknown. 

4. Results for decay constants 

We estimate the physical value of the heavy- 
light decay constants from results at (3=6.1 ob- 
tained with the 0(1 /M) NRQCD action. Since 
the value of q* is not known, we take the static 
result g*=2.18/a § as a guide, and calculate 
the central value from an average of results for 
q*=ir/a and 1/a. We then find that 

f b = 162(7)(5)(5)(11)(6)(^ 1 ) MeV, (4) 
f Bs = 190(5)(5)(5)(13)(6)(1 3 9 9 )(^) MeV. (5) 

The first error is statistical including that from 
chiral extrapolation. Remaining are systematic 
errors arising from (i) the uncertainty of q* esti- 
mated by dispersion of results for q*=ir/a and 
1/a, (ii) 0(1/ M 2 ) corrections estimated from 
comparison of results with the 0(1/ M) and the 



0(1/M 2 ) calculations, (iii) 0(a s /(aM) 2 ) errors 
estimated by dividing 0(a s /(aM)) contribution, 
which is derived from the result with static per- 
turbative correction, by aM, (iv) scaling viola- 
tion from comparison of the values at /3=6.1 with 
those at /3=5.9 and 5.7, and (v) uncertainty in 
a -1 , where the upper and lower errors correspond 
to the choice a _1 =2.62 GeV from charmonium ls- 
lp splitting and 2.21 GeV from Jk as quoted in 
Ref. Q , respectively. For fg, the central value is 
obtained with k s for strange quark fixed by mjf, 
and the last error is estimated from the shift when 
k s is derived from m^. An 0(a s A QCD /M) error 
coming from the action ([!]) is not included. A 
naive estimate of this error gives ~2% at /3=6.1. 
Some systematic errors cancel in the ratio 



SbJSb = 1.18(3)(5)(±g) 



(6) 



where the statistical error, scaling violation, and 
the uncertainty of K s , which remain, are given in 
this order. 
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